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ABSTRACT 



Context. Massive hot stars produce dense ultraviolet (UV) photon fields in their surroundings. If a very high-energy (VHE) gamma- 
ray emitter is located close to the star, then gamma-rays are absorbed in the stellar photon field, creating secondary (electron-positron) 
pairs. 

Aims. We study the broadband emission of these secondary pairs in the stellar photon and magnetic fields. 

Methods. Under certain assumptions on the stellar wind and the magnetic field in the surroundings of a massive hot star, we calculate 
the steady state energy distribution of secondary pairs created in the system and its radiation from radio to gamma-rays. 
Results. Under the ambient magnetic field, possibly high enough to suppress electromagnetic (EM) cascading, the energy of sec- 
ondary pairs is radiated via synchrotron and single IC scattering producing radio-to-gamma-ray radiation. The synchrotron spectral 
energy distribution (SED) is hard, peaks around X-ray energies, and becomes softer. The IC SED is hard as well and peaks around 
10 GeV, becoming also softer at higher energies due to synchrotron loss dominance. 

Conclusions. The radio emission from secondary pairs is moderate and detectable as a point-like and/or extended source. In X-rays, 
the secondary pair synchrotron component may be dominant. At energies S 10 GeV, the secondary pair IC radiation may be dominant 
over the primary gamma-ray emission, and possibly detectable by the next generation of instruments. 

Key words. Gamma-rays: theory - X-rays: binaries - stars: individual: LS 5039 - Radiation mechanisms: non-thermal 



1. Introduction 

Recently, several compact TeV emitters harboring an OB type 
star have been found in our Galaxy: PSR B 1259-63 is a mas- 
sive binary system containing a young non-accreting pulsar de- 
tected by HESS (Aharonian et al. 2005a); LS 5039 is likely a 
, high-mass microquasar also detected by HESS (Aharonian et al. 
■ I2005bl ): LS I -1-61 303 is a high-mass X-ray binary detected first 
; by MAGIC (Albert et al. "2006^ and recently also by VERITAS 
, (Maier 12007); Cygnus X-1 is a high-mass microquasar harbor- 
ing a black-hole detected by MAGIC (Albert et al. 120071 ) . In 
this type of objects, a significant fraction of the energy radi- 
ated above 100 GeV can be absorbed via photon-photon inter- 
actions (e.g. Ford 1984; Protheroe & Stanev 1987; Moskalenko 
& Kar akula 1994 Bednarek 199?; Boet tcher & Dermer 20()5l 
Dubus I2006al Khangulyan et al. l2007l Reynoso et al. [20081 ) 
producing secondary pairs that, under the strong radiation and 
magnetic fields present close to the massive star, can radiate ef- 
ficiently via IC or synchrotron processes. In the case in which 
IC is the dominant cooling channel, EM cascades can develop 
(see, e.g., Aharonian et al. |2006ai Bednarek 120061 Khangulyan 
et al. 120071 Orellana et al. 120071 ). If synchrotron emission is the 
dominant cooling channel, the absorbed radiation is reemitted 
at lower energies (e.g. Khangulyan et al. 120071 ) . We note that 
the scenario described here could take place not only in micro- 
quasars or non-accreting pulsar systems but in any source har- 
boring a gamma-ray emitter embedded in a dense photon field. 
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In this work, we study the radiation of the secondary pairs 
created by photon-photon absorption of gamma-rays in the sur- 
roundings of a hot massive star. We focus mainly on the case 
when the ambient magnetic field is large enough to suppress 
effectively EM cascading. The structure of the paper goes as 
follows: in Sect. |2] the adopted model is described; in Sect. [3] 
the results are presented and discussed, and our conclusions are 
given in Sect.|4] 

2. A model for pair creation and secondary pair 
emission 

2.1. The general picture 

We describe here the physical system formed by an OB star and 
a gamma-ray emitter separated by a short distance. The star is 
approximated as a point-like source of (monoenergetic) photons 
of energy eo and luminosity L,. The emitter, located at distance 
R - d» from the star and assumed to be point-like, isotropi- 
cally produces gamma-rays with an injection power Ly inj, fol- 
lowing a power-law distribution of index F. We consider here 
only the primary radiation above the minimum threshold energy 
for secondary pair production, i.e. th = l/^o (in '«eC^ units), 
since we focus on the study of the production of secondary pairs 
and their emission. The maximum energy of the primary pho- 
tons is fixed at 100 TeV. Throughout most of this work, we 
characterize the star/VHE emitter as follows: L» = 10-''' erg/s, 
dt - 3 X 10'^ cm s"', and eo = 2 x 10"^ (corresponding to a star 
temperature ^ 40000 k); we take two values for F, 2 and 3, to 
account either for hard or soft primary gamma-ray spectra. 
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Fig. 1. Sketch of the physical scenario considered here. The 
long-dashed lines represent the wind, the straight dashed lines 
the secondary pair radiation, the straight solid lines the primary 
gamma-rays, and the curved sohd hnes the advected/diffusing 
secondary pairs. 



Important ingredients of the model are the stellar wind and 
the ambient magnetic field. OB stars present fast winds with 
mass loss rates ~ 10"^ - 10"^ M© yr'^ and velocities 
y„ ~ (1 - 3) X 10^ cm s-^ (e.g. Puis et al. |2006| . We adopt 
here - 10"'' Mq yr"' and y„ = 2 x 10^ cm s"'. Moreover, 
the star generates a magnetic field in its surroundings with values 
in the stellar surface that might be as high as Bq ~ 1000 G (e.g. 
Donati et al. 2002; Hubrig et al. |2007l l. We assume here that a 
significant fraction of the magnetic field is disordered. A sketch 
of our scenario is presented in Fig.[T] 

2.2. secondary pairs in the system 

To compute the secondary pair injection spectra in the different 
regions of the binary system, we use the anisotropic differen- 
tial pair production cross section given by eq. (15) in Bottcher 
& Schlickeiser ( 119971 ). The secondary pair injection spectrum 
presents a slope similar to that of the primary gamma-rays, 
with the low-energy cutoff similar to the pair creation one, 
which depends on the angle between the two incoming photons: 
~ eth = 2/[eo(l - cos 6yy)]. The point-like approximation for the 
star fails for R ^ 2R, ~ 20 Rq from the stellar center For our 
choice of values of d^, R^, L», and eo, only photons with en- 
ergy > few TeV emitted at an angle ^ 30° with respect to the 
star-emitter line can reach the region where the finite size of the 
star comes into play (Dubus l2006al l. The monoenergetic photon 
target approximation is roughly valid provided we deal with a 
broad distribution of primary gamma-rays. The deviation from 
the real case is small and will be neglected here. 

We assume that diffusion takes place in the Bohm regime. 
Once secondary pairs are injected in different parts of the sys- 
tem, they isotropize and suffer wind advection, which is the 
dominant transport mechanism in the system. Ionization losses 
in the wind, synchrotron emission in the ambient magnetic field, 
and IC scattering of the stellar photons, are the relevant cooling 
processes of secondary pairs (y). Relativistic Bremsstrahlung of 
the secondary pairs in the stellar wind is not considered since 
the timescales for this process are much longer than ionization, 
other radiative, and advective timescales. All dependencies on 
R of the relevant physical quantities (the wind density, and the 



energy densities of the ambient magnetic and stellar radiation 
fields) are taken as oc l/R^. We consider that the TeV emitter pro- 
duces radiation long enough for the formation of a steady state 
distribution of secondary pairs in the system («(y, R)), which can 
be obtained from the following differential equation: 



dn(y,R) dj(j,RMy,R) 

Vw — — + = q{y,R), 



dR 



dy 



(1) 



where y is the particle Lorentz factor and q{y, R) is the secondary 
pair injection rate as a function of energy and distance (derived 
in the Appendix). 

A restriction on n{y,R) is that EM cascades above 
~ 100 GeV must be effectively suppressed by the ambient B, 
which must therefore be above a certain critical value, Be (see 
Khangulyan et al. |2007l l: 



lO^^erg/s 



1/2 



^1 G 

R* 



(2) 



Eq. (|2]l, together with the /^-dependence assumed here, shows 
that EM cascades are already suppressed for a w 20 G magnetic 
field atR-d,. This value appears quite moderate when looking 
at the possible stellar surface magnetic field values given above. 
However, even in the case of weaker magnetic fields, n{y,R) 
does not change much if EM cascades are not accounted for 
F ^ 2.5, since the amount of EM cascade reprocessed energy 
will be relatively small. 

For the adopted parameter values in our model, most of the 
energy of the primary VHE radiation is absorbed for distances 
^ 10'^ cm from the emitter and ~ d» from the star. Depending 
on r, radiative cooling leads to an electron distribution produc- 
ing synchrotron emission that peaks at either sharply X-ray en- 
ergies or smoothly in the range X- to soft gamma-rays. Along 
with this a fraction of energy is radiated in the GeV range via IC 
scattering. The fact that the minimum energy of the injected sec- 
ondary pairs is ~ e^in th^eC^ implies that, even under strong syn- 
chrotron or Thomson IC cooling, the emission of particles with 
energies below this value will not dominate the total radiative 
output. Ionization losses and wind advection lead to steady state 
secondary pairs radiating very little radio emission for R ^ d, 
from the star However, the wind transports secondary pair en- 
ergy to larger R, where this energy can be still efficiently ra- 
diated. In these farther regions, synchrotron radio emission can 
eventually become the dominant radiative channel. 

3. Results 

3.1. Injection and evolution of secondary pairs 

In this section we discuss the spectrum of the secondary pairs in- 
jected in the system, Qintiy) {- J qij, R)dR), and the secondary 
pair energy distribution once the steady state is reached, Mnt(T) 
{— J niy,R)dR). For illustrative purposes, we perform our cal- 



7 inj 



3 X 10^^ erg s"'. 



which is similar to the 



culations using L 
value inferred in Khangulyan et al. (120071) for LS 5039. We re- 
mark that the secondary pair emission luminosity scales linearly 
with Ly inj . 

In Figs.|2]and[3] upper panels, we show Qi^ 7^ for the whole 
volume, adopting F = 2 and 3, respectively. In the same figures, 
lower panels, A^int y^ for the whole region is also shown. In all 
the plots, the contributions to the total Qim 7^ and A^int y^ from 
R < dt and > dt are presented as well. 

As seen in the figures, Qinti7) has a similar shape to that of 
the primary gamma-rays oc er;^. As mentioned in Sect. |2] the 
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Fig. 2. Upper panel: Total QmiJ^ (solid line). Qmij^ ioxR < d» 
(long-dashed line) and > dt (dotted line) are also shown. We 
adopt r = 2. Lower panel: Total A^int (solid line). A^int for 
R < d* (long-dashed line) and > d* (dotted line) are also shown. 
We adopt F = 2, and Bo = 1 (thin lines) and 100 G (thick 
lines). The dominant cooling mechanisms at different energies 
are given. 



minimum energy of the injected secondary pairs is the pair cre- 
ation threshold energy, i.e. eth^WeC^, which changes with the an- 
gle between the two incoming photons. As a consequence, both 
qiy, R) and n(y, R) depend strongly on R as well as on P. In addi- 
tion, the shape of A^int(T) is related to the dominance of different 
cooling mechanisms and the impact of wind advection. 

For particle energies > e^in th^eC^, depending on y, Nimiy) 
is: cooled by synchrotron losses (oc y for high B in both 
R < dt and > c/»); uncooled (injection spectrum slope; for low 
B and R < dt\ because of the fast advection of the particles 
from this relatively small region); or cooled by KN IC (oc y"*^^' ; 
mainly for low B and R > dt', because there particles have time 
to lose energy). For particle energies < e,nin th'^eC^, Mnt(T) is: 
cooled by synchrotron/Thomson IC (oc y^^; R < d,); advection 
dominated (oc -y^^; R > d^; because of advection of cooled par- 
ticles created atR< d,); or cooled by ionization losses (oc con- 
stant; in both R < dt and > c/,). In the range y ~ 10 - 100, the 
1 /7?^-dependence of losses leads to a strong pileup of particles 
at/? » J,, giving a spike feature in the particle energy distribu- 
tion. The cooling mechanisms in different y ranges for F = 2 are 
pointed to in Fig. |2] lower panel; for F = 3, despite the softer 
slope,the pattern is the same. 

In Fig. m we present the total energy (Eim in Fig. |4]i accu- 
mulated in secondary pairs up to distance R, and the total energy 
rate {dEint/dt in Fig. HJi from pair creation accumulated up to R. 
The dEiaildt curve is shown up to the distance at which injection 
becomes negligeable, i.e. /?inj ~ AU. As seen in Fig.H] a signif- 
icant amount of energy, ~ lO"'^ erg, is transported up to several 
AU. Despite this energy budget, the efficiency of the emission 
from large R will be limited by the low magnetic and radiation 
energy densities and the ionization loss dominance. 
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Fig. 3. The same as in Fig. |2]but for F = 3. 
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Fig. 4. The /^-distribution of accumulated energy in secondary 
pairs (long-dashed line), and the /^-distribution of accumulated 
energy rate via pair creation (dotted line), are shown. The loca- 
tion of the emitter is marked by a vertical thick dot-dashed line, 
at \og(dt). In the plot, the star location would be at the left. 



3.2. Emission from secondary pairs 

The calculated SEDs of the synchrotron and IC emission pro- 
duced by the injected secondary pairs are shown in Figs. |5] and 
|6] In Fig. |5] we have adopted a B - 100 G, and F = 2 (upper 
panel) and 3 (lower panel). The SEDs of the emission gener- 
ated at R < dt and > dt from the star are also shown. In Fig.|6] 
the SEDs have been computed adopting B = 1 (dotted line), 
10 (long-dashed line) and 100 G (solid line), and F = 2 (upper 
panel) and 3 (lower panel). 

For energies of photons produced by secondary pairs with 
energies > emm th'WeC^, the spectrum becomes softer for larger B 
due to the dominating synchrotron channel. For B = 1 G and 
F = 2, there is a clear hardening in the synchrotron and IC spec- 
tra produced by KN IC cooling (although a proper treatment re- 
quires EM cascading). In this photon energy range, the larger F 
is, the softer the synchrotron and IC spectra are. 

For energies of photons produced by secondary pairs with 
energies < emin th^WeC^, the emission is produced either by syn- 
chrotron/Thomson IC cooled secondary pairs (R ^ dt), or by 
secondary pairs affected already by wind advection {R > dt). 
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The spectral hardening at low energies in the synchrotron and 
IC emission for R ^ is produced by ionization losses. The 
particle pileup around y ~ 10 impacts on the low energy radia- 
tion from R> d». 

The inhomogeneity of q{y, R), and the sensitivity of n(y, R) 
to the cooling and transport conditions, imply that the produced 
emission changes strongly with location in the system. Further 
complexities like a magnetic field with angular dependences, 
or a strongly inhomogeneous diffusion coefficient, would yield 
globally and locally different SEDs. 

3.2.1 . LS 5039 and Cygnus X-1 

Adopting Z?o ~ 100 G, our model predicts radio fluxes and X- 
ray luminosities of several mJy and several 10^^ erg s"', respec- 
tively. The parameter values adopted in this work are similar to 
those of LS 5039 and Cygnus X-1, making worthy a discussion 
of the results in the context of these two sources. 

The obtained radio fluxes and X-ray luminosities are not far 
from those found in LS 5039 (e.g. Paredes et al. 2000; Bosch- 
Ramon et al. 200 7 ), implying that the secondary pair contribu- 
tion in this source may be comparable, if not dominant, to that of 
any intrinsic radio and X-ray component linked to the TeV emit- 
ter itself. We remark that, especially in the radio band, a slightly 
higher B-field would yield significantly larger radio fluxes due 
to the B^-dependence of the synchrotron emission. The moderate 
X-ray luminosity in LS 5039 permits a look for such a secondary 
pair component. In the GeV regime, the secondary pair IC fluxes 
are far from those found by EGRET (Paredes et al. 2000 ), sug- 
gesting that the latter are due to the combination of both primary 
and secondary pair IC radiation. 

In the case of Cygnus X-1, adopting the same Lyinj as for 
LS 5039, the predicted radio flux from secondary pairs is a sub- 
stantial part of the total radio flux of the source (e.g. Stirling et al. 
1200 11 1. In X-rays, explained as thermal emission from accre- 
tion disk/corona-like regions, the high luminosity of Cygnus X-1 
makes it impossible observationally to disentangle a possible 
secondary pair component. 

The lack of evidence of accretion in LS 5039 (Bosch-Ramon 
et al. |2007l and references therein) is an important difference be- 
tween this source and Cygnus X-1, which shows clear X-ray 
accretion features (Albert et al. 12007! and references therein). 
Nevertheless, as noted, e.g. by Bogovalov & Kelner ( |2005l l, the 
fact that in some sources the accretion disk luminosity is unde- 
tectable does not imply lack of accretion. This may explain why 
LS 5039 does not show accretion signatures in its X-ray spec- 
trum. 

3.2.2. LS I +61 303 

For a compact system with the properties listed in Sect. |2] only 
VLBI interferometers with angular resolution $0.1 milliarcsec- 
onds would be able to resolve the radio emission. Nonetheless, 
within the constraints of our model, we have explored the pos- 
sibility of explaining the extended radio emission of a less 
compact system, LS I +61 303, which presents a peculiar ex- 
tended radio structure with changing morphology along the orbit 
(Dhawan et al. |2006l l. We have computed the radio emission pro- 
duced at the spatial scales of the observed extended radio struc- 
tures. Adopting B{) - 100 G, we have obtained the SEDs for 
the synchrotron radio emission originated in different regions: 



' The relevant parameters characterizing the system properties can be 
found in Bosch-Ramon et al. ( 12006b . 




log ( E leVJ ) 

Fig. 5. Upper panel: Computed SEDs for the synchrotron and 
the IC secondary pair emission. We adopt Z?o = 100 G and F = 
2. The emission from the whole volume (solid line) filled by 
particles, and that of a region located at R < d, - ~ 0.1 mas at 
3 kpc- (long-dashed line) and > dt (dotted line). Lower panel: 
The same as upper-panel but for F = 3. 



R < 1, > 1, > 2, and > 3 AU. The radio and broadband (syn- 
chrotron plus IC) SEDs are presented in Fig.|7](lower -left- and 
upper panels, respectively); the spatial distribution of the cor- 
responding emitting particles around LS I +61 303 after an in- 
jection time of one orbital period is also shown (lower -right- 
panel). 

To implement these calculations, we adopt the orbital dis- 
tance corresponding to the phases when the source was detected 
by MAGIC (Albert et al. |2006] l: ~ 6 x 10'^ cm. F is taken 2.6 
(Albert et al. I2006I I. and Ly i„j = 3 x 10''-'' erg s"', the same 
as the one taken for LS 5039, enough to explain observations. 
We note that LS I +61 303 is not detected by MAGIC above 
a few hundred GeV when outside the phase range ~ 0.5 - 0.7. 
Nevertheless, primary gamma-rays with softer spectra may still 
be injecting a significant amount of secondary pairs in the sys- 
tem despite being barely detectable above few hundreds GeV. 
For the spatial distribution of particles, we followed individu- 
ally their energy and spatial evolution accounting for advection 
and (Bohm) diffusion in the wind using a simple Monte-Carlo 
simulation (Bosch-Ramon et al., in preparation). The orbital pa- 
rameters of LS I +61 303, like the eccentricity of the system 
(e = 0.72; see Casares et al. ,2005, for the system parameters), 
were considered. 

From Fig.|7] it is seen that radiation flux levels of ~ 20 mJy 
(8 GHz) are reached. In this source, most of the secondary pair 
radio emission would appear extended, pointing in the direc- 
tion opposite to the star and bending due to orbital motion, as 
shown in the secondary pair spatial distribution in Fig. |7] This 
kind of behavior is similar to that found by Dhawan et al. (I2006I I 
in LS 1 61 +303. These authors associated the radio morphology 
of this source to a particular non-accreting pulsar scenario (see 
Dubus |2006bl l, although recent hydrodynamical simulations of 
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Fig. 6. Upper panel: Computed SEDs for the synchrotron and 
the IC secondary pair emission. We adopt Bq - I (dotted line), 
10 (long-dashed line) and 100 G (solid line), and F - 2. Lower 
panel: The same as upper panel but for F = 3. 



stellar/pulsar colliding winds predict quite different morpholo- 
gies (see Romero et al. 2007 and Bogovalov et al. 120071) . 

The predicted X-ray luminosities are similar to those ob- 
served in this source (e.g. Sidoli et al. I2006i) . which, as in 
LS 5039, is small enough to allow the study of a secondary pair 
component. Despite the fact that the complexity of the X-ray 
emission cannot be completely explained by our model, the pro- 
cesses considered here may be behind a significant fraction of 
the observed X-ray emission. Like LS 5039, LS I H-61 303 was 
also proposed to be a (variable) GeV source (Tavani et al. 1 19981 ). 
and the secondary pair contribution to this energy range may be 
significant (see also Bednarek ,2006i ). 

4. Conclusions 

We conclude that the presence of a powerful VHE emitter near 
a massive hot star leads unavoidably to non-thermal emission in 
the stellar wind. The question whether this secondary pair radi- 
ation is detectable depends on the magnetic and radiation field 
strength in the system. 

In the scenario explored here, we predict moderately hard 
spectra, and fluxes of several mJy for typical galactic distances, 
not far from the fluxes detected in some microquasars. We re- 
mark that radio emission is produced mainly in the regions with 
R > d» and may be dominant over any other radio component 
(e.g. linked to the TeV emitter itself), and the radio morphology, 
flux and spectrum are strongly sensitive to the geometry and R- 
dependence of the magnetic field. Interestingly, the predicted ra- 
dio morphology could resemble that shown by the VHE emitting 
X-raybinai-yLSI-H61 303. 

The fluxes of the (likely non-thermal) X-ray emission of sev- 
eral TeV emitting X-ray binaries (e.g. LS I -h61 303, LS 5039), 
typically around lO^^"' - lO"''* erg s"', are similar to the values 
predicted here. It implies that the secondary pair X-ray radia- 
tion could be comparable to, or even dominate over, any intrin- 



sic component linked to the TeV emitter itself. It is worth not- 
ing that due to their moderate X-ray luminosities LS 5039 and 
LS I +61 303 are good candidates to look for a secondary pair 
contribution in this energy range. 

We show here that the secondary pair emission at ~ GeV 
energies, if not strongly dominated by an intrinsic GeV emit- 
ter, could be revealed by GLAST. Above e^in th ~ 10 GeV, soft 
primary spectra and/or moderate-to-high magnetic fields would 
imply low fluxes. In the latter case, significant synchrotron en- 
ergy losses would suppress EM cascading. 

Our calculations show that even with simple assumptions 
on the system geometry, primary gamma-ray injected spectra, 
and B-field and wind structure, the resulting situation is quite 
complex. Thus a detailed characterization of the secondary pair 
non-thermal emission in a particular source is a difficult task for 
which stellar wind physics, primary VHE emission modeling, 
and high quality data are required. 
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Fig. 7. Upper panel: Computed SEDs for the synchrotron and the IC secondary pair emission produced in the whole volume (solid 
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of the emitting particles, for the orbital plane (a) and a plane perpendicular to it (b), in the surroundings of LS I +61 303 after one 
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5. Appendix 

The secondary pair injection function for a monoenergetic point- 
like source of target photons is determined by the following in- 
tegral: 



1 



dcos^ ■ 



(l-cosQ'(r)) dA^^_^do^ 
dedf dy 



de- 



(3) 



where dA'y(e)/dedf is the primary gamma-ray injection spectrum 
per time unit, r/r is the distance/optical depth from the gamma- 
ray emitter to the secondary pair creation location, and a{r) is 
the interaction angle at this location. A sketch of the situation 
is presented in Fig. [8] The cross-section do-p/dy is given be the 
eq. (15) from Bottcher & Schlickeiser (1997) and the kinematics 
constraints: 



<y < - 



+ ef„(l-cosa-(r))) 



(4) 



2 ' 2 

The distance and the interaction angle are defined as the follow 
ing: 

+ - dl 



■ R^ +dl 



2d,R cos , cos a{x) 



2px 



(5) 
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Fig. 8. Geometry of the photon-photon interaction. 



where p = -ijx^ + dl + x{R^ - dl - r^)/r is the distance from 
the optical star to the gamma-ray absorption point. Finally, the 
optical depth t is calculated like 



T(r,R) 



r 

I 



(1 - COS a(x)) 
dx 2 o-p ' 



(6) 



where CTp is the total pair production cross-section (see e.g. 
Berestetskii et al. ( |1982| i). 



